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Abstract 

In  this  paper  a  fractal  permeability  model  for  the  gas  diffusion  layer  (GDL)  of  PEM  fuel  cells  (PEMFCs)  is  presented.  The  model  accounts  for 
the  actual  microstructures  of  the  GDL  in  terms  of  two  fractal  dimensions,  one  relating  the  size  of  the  capillary  flow  pathways  to  their  population 
and  the  other  describing  the  tortuosity  of  the  capillary  pathways.  In  addition,  the  gas  molecule  effect  is  considered  by  using  the  Adzumi  equation. 
The  fractal  permeability  model  is  found  to  be  a  function  of  the  tortuosity  fractal  dimension,  pore  area  fractal  dimension,  sizes  of  pore  and  the 
effective  porosity  of  porous  medium  without  any  empirical  constants,  mercury-intrusion  porosimetry  was  used  to  measure  the  microstructures  of 
the  GDL.  Based  on  scanning  electron  microscope  (SEM)  images,  two  fractal  dimensions  are  determined  by  the  box-counting  method.  To  verify 
the  validity  of  the  model,  the  predicted  permeability  data  of  the  present  fractal  model  were  compared  with  the  experimental  data  supplied  by  Toray 
Inc.  It  is  found  that  the  permeability  prediction  of  the  model  was  in  accordance  with  experimental  data.  This  verifies  the  validity  of  the  present 
fractal  permeability  model  for  the  GDL. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

PEM  fuel  cells  (PEMFCs)  have  a  high  efficiency,  power 
density,  reliability  and  a  quick  startup  capability.  They  could 
have  a  wider  use  in  spaceflight  and  military  applications  with  a 
huge  market  potential  for  fuel  cells  electric  utilities,  automotive 
transportation  and  in  portable  equipment.  In  the  simulations  of 
PEMFCs,  supposing  that  porous  media  including  gas  diffusion 
layer  (GDL),  catalyst  layer  and  membrane  are  always  homoge¬ 
neous,  diffusion,  momentum  and  energy  equations  are  solved 
together  with  macroscopic  transport  properties  such  as  perme¬ 
ability,  effective  thermal,  electronic  conductivity  and  effective 
diffusion  coefficient.  Firstly,  this  assumption  is  wrong.  In  real¬ 
ity,  the  GDL,  catalyst  layer  and  the  membrane  are  anisotropic. 
Secondly,  it  relates  to  the  transport  mechanisms  of  PEMFCs 
without  consideration  of  the  micro-geometry  and  connectivity 
of  the  chaotic  pore  spaces  of  actual  porous  media  through  which 
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the  transport  takes  place,  and  such  simulations  give  limited  direc¬ 
tion  for  the  correct  design  of  PEMFCs.  Since  porous  media  are 
complex,  there  is  no  good  way  to  represent  their  microstructures 
exactly.  Thompson  et  al.  [  1]  confirmed  that  the  interspaces  in  real 
porous  media  have  a  fractal  character  through  scanning  electron 
microscope  (SEM).  This  work  tries  to  apply  fractal  theory  to 
PEMFC  modeling,  and  develops  a  fractal  model  of  permeabil¬ 
ity,  one  of  the  transport  characteristics  of  the  GDL.  This  model 
considers  the  microstructures  of  porous  media,  and  will  help  us 
better  understand  the  heat  and  mass  transport  mechanisms  in 
this  porous  layer  of  the  PEMFC. 

2.  Research  methods  and  fractal  research  of 
permeability 

In  this  section,  conventional  research  methods  and  fractal 
research  progress  on  permeability  are  introduced. 

Research  methods  of  the  permeability  included  experimen¬ 
tal  measurement,  analytical  solution  and  numerical  simulation. 
Traditionally,  determination  of  the  permeability  tensor  is  best 
accomplished  through  experimental  measurements  since  many 
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and  various  measurement  techniques  all  yield  consistent  results. 
However,  these  measurements  often  require  a  large  number  of 
carefully  controlled  experiments  and,  in  general,  have  no  pre¬ 
dictive  capability  (i.e.,  each  new  material  must  be  handled  on 
a  case-by-case  basis).  As  a  consequence,  analytical  (for  simple 
cases)  or  numerical  predictions  of  the  permeability  have  for  the 
past  few  years  been  the  goal  of  many  researchers  in  the  field  of 
porous  media  flow. 

The  permeability  prediction  models  are  mainly  of  three  types: 


(a)  Fitting  the  experimental  data  to  a  pure  empirical  formula 
with  several  empirical  constants.  These  empirical  constants 
usually  do  not  have  any  specific  physical  meaning,  and  the 
constants  gained  by  different  people  always  differ  in  values, 
therefore  lacking  of  consistency.  For  example: 

Lei  et  al.  [2]  studied  the  relation  between  porosity  and 
permeability  of  the  grain  bed,  the  empirical  formula  gotten 
from  the  experiment  data  is 

K  =  1.72  x  l04dlA65<p4'69  (1) 


where  d  and  0  are  the  grain  diameter  and  porosity,  respec¬ 
tively,  Eq.  (1)  is  applied  to  the  range  of  d=  0.1-0.45  mm. 
There  are  three  empirical  constants:  1.72  x  104,  1.465  and 
4.69,  so  it  is  obvious  that  they  have  no  specific  physical 
meaning. 

Adler  and  Thovert  [3]  presented  several  empirical  for¬ 
mulae  for  the  permeability  of  real  porous  media,  e.g.  an 
empirical  formula  for  a  glass  particle  porous  media  is 

K  =  O.Wlrlp-51  (2) 


where  re  is  the  equivalent  radius  of  glass  particle,  Eq.  (2)  is 
applied  to  the  range  from  0  =  0.4  to  0  =  0.79.  There  are  also 
several  empirical  constants  in  Eq.  (2). 

(b)  By  specific  ideal  models  and  theoretical  analysis,  in  which 
the  relation  between  permeability  and  different  parameters 
of  the  structures  in  the  porous  media  can  be  derived,  which 
is  semi-empirical  approach.  Some  numerical  coefficients  of 
these  models  must  be  confirmed  by  experiment  for  a  specific 
porous  medium. 

The  famous  Kozeny-Carman  equation  is  a  typical  exam¬ 
ple  of  semi-empirical  formula,  given  by 


0» +  i 

C(1  -  0f 


where  the  exponential  n  and  constant  C  are  called 
Kozeny-Carman  constants,  these  two  constants  vary  for  dif¬ 
ferent  porous  media. 

(c)  Pure  theoretical  formulas  are  derived  from  Darcy’s  law. 
They  are  only  applied  to  an  ideal  structure  (such  as  circi- 
nal  capillary,  rectangle  channel,  etc.)  porous  media,  but  not 
to  the  random  and  disorderly  structure  of  real  porous  media. 


In  practical  applications,  hardly  any  porous  media  have  an 
ideal  structure.  The  transport  character  parameters  of  the  porous 
media  have  been  measured  mainly  by  experiments  or  com¬ 
puted  by  numerical  simulation  over  several  decades,  so  how 


to  decrease  or  eliminate  these  empirical  constants  is  the  subject 
that  people  in  this  field  still  work  on. 

A  large  number  of  reports  indicate  that  the  micro  structures 
and  pore  size  distributions  of  porous  media  have  a  fractal  char¬ 
acteristic.  Therefore,  it  is  possible  to  obtain  the  permeability 
of  porous  media  through  fractal  analysis  of  the  pore  structures. 
Many  people  have  made  great  efforts  to  study  fractal  theories 
and  develop  their  application;  so  far  they  have  already  made 
some  progress. 

Tyler  and  Wheatcraft  [4]  combined  porous  capillary  channel 
models  with  the  Van  Genuchten  empirical  formula,  and  got  the 
fractal  model  for  the  water  content  characteristic  curve.  Rieu  and 
Sposito  [5]  introduced  theoretical  models  of  pore  distribution 
based  on  a  Sierpinski-type  gasket.  Perfect  et  al.  [6]  improved 
models  put  forward  by  Tyler,  Wheatcraft,  and  Rieu,  Sposito, 
and  they  got  the  relation  between  the  water  content  and  the  pore 
pressure  in  porous  media.  Shepard  [7]  applied  the  Koch  curve  to 
the  simulation  of  soil’s  pore  distribution,  and  then  got  a  hydraulic 
conductivity  formula  in  the  form  of  an  exponential. 

Adler  and  Thovert  [3]  made  a  numerical  simulation  analy¬ 
sis  of  the  transport  problems  of  fractal  objects  with  geometrical 
seep  models  in  early  times.  But  the  results  did  not  have  any 
connection  with  the  fractal  dimension,  nor  were  they  compared 
with  experimental  data.  They  were  only  connected  with  geo¬ 
metrical  iterative  times.  Adler  studied  the  transport  character 
in  porous  media  later,  and  considered  that  the  permeability  of 
porous  media  could  be  expressed  as  follows:  K=K(cp ,  Df,  . . .), 
where  Df  is  the  fractal  dimension.  But  the  quantitative  expres¬ 
sion  of  permeability  was  not  given  in  this  paper. 

Chen  and  Shi  [8]  raised  a  fractal  permeability  model  of 
a  porous  soil  media,  and  this  model  has  four  empirical  con¬ 
stants.  The  area  fractal  dimension  of  the  soil  pore  distribution 
reported  in  their  paper  is  that  Df  <  1 .  But  according  to  basic  frac¬ 
tal  theories,  1  <  Df  <  2  in  two-dimensional  space  and  2  <  Df  <  3 
in  three-dimensional  space,  so  their  results  about  the  area  frac¬ 
tal  dimension  is  doubtful.  Another  shortcoming  of  their  work 
is  that  they  did  not  compare  theoretical  models  to  experimental 
results,  so  the  validity  of  the  models  needs  to  be  confirmed  fur¬ 
ther.  In  their  recent  work,  they  improved  the  computation  of  the 
fractal  dimension  and  got  the  fractal  dimension  Df  between  1 
and  2,  at  the  same  time,  the  empirical  constants  decreased  from 
4  to  2. 

Pitchumani  and  Ramakrishnan  [9]  presented  a  fractal  analyt¬ 
ical  model  for  evaluating  the  permeability  of  a  porous  fibrous 
textile.  This  model  has  not  any  empirical  constants,  which  is  a 
great  improvement,  but  their  model  is  incorrect  for  disobeying 
the  basic  fractal  theories. 

Yu  and  Lee  [10]  deduced  a  general  analytical  solution  model 
for  a  kind  of  porous  fibrous  textile 


Hg/\2  +  JkI/2 
A 


where  Hg,  A  and  k±  are  biggest  pore  area,  pore  height, 
unit  area  and  transverse  permeability  of  fasciculi,  respectively. 
Although  Eq.  (4)  has  no  empirical  constants,  it  is  not  applicable 
to  porous  media  with  random  and  irregular  structures  in  which 
flow  cannot  be  simplified  into  one-dimensional  Stokesian  flows. 
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In  their  subsequent  research  work,  a  universal  fractal  analytical 
solution  model  of  permeability  for  porous  media  was  derived 
based  on  the  Hagen-Poiseulle  equation  and  Darcy’s  law.  Their 
research  objectives  are  particle  porous  media,  such  as  grit  and 
soil,  so  their  model  does  not  apply  to  the  GDL  of  a  PEMFC  that 
is  a  fibrous  and  microporous  medium. 

In  the  following  sections,  we  will  first  introduce  the  basic 
fractal  characteristic  analytical  theories  of  the  GDL  in  a  PEMFC. 
Then  the  analytical  solution  model  of  permeability  will  be 
derived  based  on  experiment  research.  Finally  we  will  assess 
the  validity  of  the  models. 

3.  Fractal  permeability  model  for  GDL  of  PEMFC 

We  will  introduce  some  basic  concepts  of  the  fractal  theory 
first. 

Scale  relationship  between  the  measurement  of  a  fractal 
object  M(L)  and  the  metrical  yardstick  L  is  as  follows  [11]: 

M(L)  ~  LDt  (5) 

where  Df  is  the  fractal  dimension,  M{L)  the  quality,  or  volume, 
or  area  of  the  object  or  the  curvilinear  length  and  L  is  the  scale. 

Another  character  of  the  fractal  object  (such  as  pore,  or 
islands  in  the  Earth)  is  that  the  cumulative  number  N  (for  exam¬ 
ple,  the  number  of  pore)  of  the  fractal  object  and  pore  size 
distribution  are  related  by  the  equation  [11] 

N(L>X)=(^Y  (6) 

where  A  and  Amax  are  the  pore  size  and  the  maximum  pore  size, 
respectively,  as  far  as  porous  media  are  concerned.  The  first 
derivative  of  Eq.  (6)  with  respect  to  A  can  be  written  as 

-dN  =  DfX°<X-(D<+1)  d\.  (7) 


character,  so  fractal  theories  are  effective  in  the  analysis  and 
evaluation  of  the  pore  distribution  characteristic  of  the  GDL. 

In  simulation  of  a  PEMFC,  the  gas  permeability  diffusion  in 
a  GDL  needs  to  be  taken  into  account.  According  to  the  exper¬ 
imental  data  from  mercury  porosimetry,  it  can  be  seen  that  the 
pore  diameters  of  the  GDL  range  from  10-5  to  10-8  m,  and 
most  pore  diameters  are  about  10-6  m.  However,  the  mean  free 
path  of  a  gas  molecule  is  about  10-7  m;  thus,  the  Knudsen  num¬ 
ber  ranges  from  about  10-2  to  10,  so  the  slipstream  of  the  gas 
molecule  should  be  considered.  The  flow  rate  q(X)  through  a 
single  tortuous  capillary  can  be  given  by  the  Adzumi  equation 

[13] 


it  A p  A4  s  I 2ttRT  A3  A p 

128  Lx(k)  +  6  V  M  TJx)~p 


where  fi  is  the  viscosity  of  the  fluid,  A p  the  pressure  gradient 
and  s  is  the  Adzumi  constant  which  is  a  non-dimensional  pro¬ 
portional  fraction,  it  takes  a  value  of  0.9  for  one  kind  of  gas  and 
a  value  of  0.66  for  a  gaseous  mixture.  M  is  the  gas  molecular 
weight,  and  the  gas  constant  R  is  8.3143  J  mol-1  K-1. 

Eqs.  (6)— (9)  form  the  basis  of  the  present  fractal  permeability 
model  for  the  GDL  of  a  PEMFC,  which  we  will  derive  as  follows. 

The  total  flow  rate  Q  can  be  obtained  by  integrating  the  indi¬ 
vidual  flow  rate  q( A)  over  the  entire  range  of  pore  sizes  from  the 
minimum  pore  Amin  to  the  maximum  pore  Amax  in  a  unit  cell. 
According  to  Eqs.  (7)  and  (9),  we  have 


^max  rk 

Q  =  ~  I  q(X)dN(X)  =  j 


max 


7i  A p  Az 


'■min 


min 


128  Lt(A)  jl 


el  2nRTJ^Ap]  (Df+l) 

6  V  M  Lt(X)  p  1 


(10) 


The  number  of  pores  whose  sizes  fall  within  the  infinitesimal 
range  A  to  A  +  dA  is  given  by  Eq.  (7),  and  — dA>  0. 

When  liquids  flow  through  the  pores  of  porous  medium,  the 
capillaries  may  be  tortuous.  These  tortuous  capillaries  could  be 
expressed  by  fractal  equation  [12] 

Lt(X)  =  Lfx1-01  (8) 

where  Dj  is  the  tortuosity  fractal  dimension,  and  lies  in  the 
range  1  <  Dt  <  2,  which  represents  the  extent  of  convolutedness 
of  capillary  pathways  for  fluid  flow  through  a  medium.  Note  that 
for  a  straight  capillary  path  Dj=  1 ,  and  a  higher  value  of  Dt 
corresponds  to  a  highly  tortuous  capillary.  The  limiting  case  of 
Dx  =  2  corresponds  to  a  highly  tortuous  line  that  fills  a  plane.  Let 
the  diameter  of  a  capillary  in  the  medium  be  A  and  its  tortuous 
length  along  the  flow  direction  be  Lt(A).  Lo  is  representative 
length  of  channels.  For  a  straight  capillary,  Lt(A)  =  Lo. 

In  early  research,  we  made  mercury  porosimetry  and  BET 
absorption  experiments  for  the  GDL  and  from  the  experimental 
data  of  mercury-intrusion  porosimetry,  two  fractal  dimensions 
were  computed.  The  calculated  results  showed  that  linear  regres¬ 
sion  of  the  experimental  data  is  obvious  and  GDL  has  a  fractal 


Substituting  the  definitions  of  Lt(A)  into  Eq.  (10),  we  can  get 


max 


Q  = 


'-min 


it  ApDjX^  :  2+nT-nf 
128  p,L(p 


+  6 


e  [2ttRT  A p  DfXf^  l+Dl_D{ 


M  p  Lf 
n  ApDfX?*  1 


dA 


128^^0  T  3  +  Dt  —  Df 


X 


3-\-Dj— Df 


1 


k 


max 


k 


mm 


e  2nRT  A p  Df/-”[1X 

6  V  M  p  L ^  2+Dj-Df 


^2+Dj—Df 


k 


max 


k 


mm 


7Z  ApDfXPt 


max 


i 


128 3  +  Dt-  Df 


^3+Dj  —  Df 


max 


X 


1  - 


A 


A 


mm 


max 


3+DT— Df-|  / 2lzRT  Ap  D ^ 

+  6  V  M  p  Lqt 


1 


X 


_ t,  2+Dj  —  Df 

2  +  Dj  —  Df  max 


1  - 


A 


2+Dj—Df- 


min 


A 


max 


(ii) 
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Using  the  Darcy’s  laws,  the  permeability  of  GDL  could  be 
written  as 


K  = 


ijLqQ 
A  pA 


n  Ly°T 


Df 


X 


1  - 


_ y  3+Dt 

128  A  3  +  Dj  —  Df  max 

Vin  \  3+Dl~Df]  e  [2 nRT  ii  Ly°T 


X 


max 


X 


Df 


_ y  2+Dj 

2  +  Dj  —  Df  max 


1  - 


+  6 
X 


M  A  p 

2-\-Dj— Df~ 


mm 


X 


max 


(12) 


conditions.  Thus,  for  flow  through  the  unit  cell  with  a  single 
capillary  tube  or  pore  with  Df  =  2,  we  have  the  maximum  flow 
rate  and  maximum  permeability  from  Eqs.  (17)  and  (18), 


max 


ttA  p  4  _  Dg  A  A p 

128/iLo Amax  ~~  32  Loti 


max 


_2_i  4  d| 

128  A  max  32 


(19) 

(20) 


Let  the  first  item  of  Eq.  (12)  be  K\ ,  the  second  item  ,  then 
Eq.  (12)  becomes 

K  =  /G  +  K2  (13) 


In  the  two-dimensional  space,  due  to  1  <  Dt  <  2  and 
l<Df<2,  the  exponents  of  Amin/Amax,  3  +  Dt  —  Df>l  and 
2  +  Dx  —  Df  >  1,  respectively,  and  because  Amin/Amax  <  10-2,  Eq. 
(12)  could  be  reduced  to 


K  = 


pLoQ_ 

A  pA 


7T  Ly°T 


128  A 


1 3+Z>r 

3  +  Dx-Df  max 


£ 

+  6 


2nRT  /I  LpDj  Df 


M  A  p 
Thus,  K\  can  be  gained  by 

Df 


_ y  2+Dj 

2  +  Dx-Df  max 


K  i  = 


128  A 


_ y  3+Dj 

3  +  Dx-Df  max 


(14) 


(15) 


Eq.  (19)  is  exactly  the  Hagen-Poiseulle  equation  [13],  and  the 
permeability  value  of  D2/ 32  is  exactly  the  expression  for  flow 
through  a  pipe. 

K2  in  Eq.  (13)  is  the  modified  item  when  the  Knudsen  effect 
is  taken  into  account,  so  the  effect  of  K2  on  the  total  permeability 
is  not  very  large,  while  the  first  item  K\  plays  an  important  role 
in  permeability  computation.  From  Eqs.  (19)  and  (20),  it  can 
be  seen  that  our  model  is  consistent  with  the  physical  situation. 
Therefore,  one  can  find  the  flow  rate  and  the  permeability  for 
flow  through  the  unit  cell  with  a  straight  capillary  tube  either 
from  the  modified  Hagen-Poiseulle  equation  or  from  Eqs.  (19) 
and  (20)  without  consideration  of  the  Knudsen  effect. 

As  mentioned  above,  we  can  see  that  in  Eq.  (14),  K\  is  the 
permeability  expression  that  can  be  derived  from  the  modified 
Hagen-Poiseulle  equation  and  K2  is  the  modified  item  caused 
by  the  molecule  slipstream. 

Subsequently  the  expression  of  the  total  area  of  the  unit  cell 
A  will  be  derived  by  the  definition  of  porosity  as  follows: 


According  to  Eq.  (11),  the  corresponding  total  flow  rate  Q\ 
can  be  written  as 


X  1  1  3+Dt 

128 /iLjf  3  +  DT  -  Df  max 


(16) 


Eq.  (15)  indicates  that  the  permeability  K\  is  a  function  of 
the  pore  area  fractal  dimension  Df ,  the  tortuosity  fractal  dimen¬ 
sion  Dx  and  structural  parameters  A,  Lo  and  Amax.  If  a  straight 
capillary  model  (Dx  =  1)  is  assumed,  Eqs.  (15)  and  (16)  can  be 
reduced  to 


it  A p  Df 


A 


4 


128/xLo  4  —  Df  max 
1  Df 


71 


128  A  4  -  Df 
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max 


(17) 

(18) 


respectively.  Eqs.  (17)  and  (18)  indicate  that  the  flow  rate  Q\ 
and  permeability  K\  are  very  sensitive  to  the  maximum  pore 
size  Amax.  It  is  also  shown  that  a  higher  value  of  Df  corresponds 
to  higher  flow  rate  and  higher  permeability  values.  From  Eqs. 
(17)  and  (18),  we  can  see  that  the  flow  rate  and  the  permeability 
will  reach  the  possible  maximum  values  as  the  pore  area  fractal 
dimension  approaches  its  possible  maximum  value  of  2.  The 
limiting  case  of  Df  =  2  corresponds  to  a  smooth  surface  or  a 
plane  or  a  compact  cluster.  This  means  that  if  we  consider  the 
smooth  surface,  compact  cluster,  circle  or  square  to  be  the  cross- 
section  of  a  pore,  the  fractal  dimension  of  the  cross-section  is 
2  and  the  pore  volume  fraction  of  the  cross-section  is  1.  Both 
the  flow  rate  and  the  permeability  are  a  maximum  under  such 
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where  Ap  is  the  pore  area  of  the  unit  cell. 

According  to  the  definition  of  porosity,  we  get 


71  ^f^  max 
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^  A  0  40(2 -Df) 


1  - 


X 


2-Dn 


mm 


X 


max 


(22) 


where  0  is  the  porosity. 

Substituting  Eq.  (22)  into  Eq.  (12),  we  have 
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In  the  two-dimensional  space,  Eq.  (23)  can  be  abbreviated  to 
_  i  2-pf 

32  3  +  DX  -  Df  [1  -  (VinAmax)2“Df] 

e  j2RTixLf)f  (4-2 D{)4> _ kf*, 

~3  V  JtM  p  2  +  Dt  -  Df  [1  -  (WW)2_Df] 

(24) 


Eqs.  (23)  and  (24)  without  empirical  constants  are  the  ana¬ 
lytical  solution  model  of  permeability  for  GDL.  It  can  be  seen 
that  the  permeability  K  is  a  function  of  two  fractal  dimensions 
Df,  Dt  and  porous  material  microstructure  parameters.  Eqs. 
(23)  and  (24)  all  have  1 1  parameters,  and  every  parameter  has 
specific  physical  meaning.  So  the  permeability  model  gained 
by  the  fractal  theories  is  obviously  better  than  the  well-known 
Konzeny-Carman  equation  from  the  viewpoint  of  describing  a 
penetrative  mechanism. 

According  to  the  following  formula  proved  recently  by  Yu 
and  Li  [14]: 
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mm 


A 


max 


dE—Df 


(25) 


where  d e  is  Euclid  dimension,  and  =  2,  1  <  Df  <  2  in  the  two 
dimensional  space;  while  d^  =  3,  2  <  Df  <  3  in  the  three  dimen¬ 
sion  space.  Eq.  (24)  can  be  further  abbreviated  to 

K=J  2 -Of 

32  3  +  Dt  —  Df  1  -<p 


s  /2RT^Lfff  (4  —  2Df)0  k^ 
+  3  M  ttM  p  2  + Dj  -  D{\  -  & 


(26) 


In  the  following  section,  we  will  discuss  the  determination 
of  two  fractal  dimensions  Df  and  Dt  and  fractal  permeability 
prediction  of  GDL. 


4.  Two  fractal  dimensions  evaluation  and  permeability 
prediction  for  the  GDL 


Table  1 


The  microstructual  parameters  and  operating  parameters  of  samples 


Parameter 

Number  one 

Number  two 

Description 

An  ax 

8  x  10-5  m 

7  x  10-5  m 

Maximum  pore 
diameter 

kmin 

3.079  x  10"8  m 

1.487  x  10“8m 

Minimum  pore 
diameter 

0 

0.55 

0.78 

Porosity 

Lo 

1.9  x  10-4  m 

1.9  x  10-4  m 

Gas  diffusion  layer 
thickness 

li 

1.2  x  10-5  Pas 

1.2  x  10-5  Pas 

Hydrogen  viscosity 

P 

1.01325  x  105  Pa 

1.01325  x  105  Pa 

Hydrogen  pressure 

T 

353  K 

353  K 

Temperature 

samples  Dt  =  1.1447,  the  pore  area  dimension  of  TGP-H-060 
carbon  paper  Df  =  1.9669  and  that  of  TGP-H-060  carbon  paper 
treated  with  PTFE  Df  =  1.9276. 

As  follows,  the  prediction  of  permeability  will  be  made  when 
H2  goes  through  two  samples  (number  one  sample  is  TGP-H- 
060  carbon  paper  and  number  two  sample  is  TGP-H-060  carbon 
paper  treated  with  PTFE).  Table  1  gives  the  operating  parameters 
and  some  experimental  data  from  mercury  porosimetry. 

We  now  compare  the  permeability  values  based  on  the  present 
fractal  model  with  the  reference  value  supplied  by  Toray  Inc.  for 
the  number  one  sample.  The  results  are  presented  in  Fig.  1. 

The  real  line  in  Fig.  1  shows  that  the  fractal  permeability 
model  for  the  number  one  sample  is  slightly  lower  than  the  refer¬ 
ence  value,  that  is,  8  x  10~8  m2.  In  the  same  figure,  we  also  plot 
the  well-known  Kozeny-Carman  equation  given  by  Eq.  (3),  in 
which  C  =  180/A^ean,  0  =  0.78,  the  mean  pore  diameter  Amean 
is  computed  by  experimental  data  from  mercury  porosimetry.  It 
is  shown  that  the  present  fractal  permeability  model  is  in  better 
agreement  with  the  reference  value  than  the  Kozeny-Carman 
equation.  Moreover,  we  forecasted  the  permeability  of  TGP-H- 
060  carbon  paper  treated  with  PTFE  by  our  fractal  model  and 
Kozeny-Carman  equation,  respectively,  also  shown  in  Fig.  1. 


In  the  fractal  analytical  solution  model  of  permeability,  we 
need  to  determine  two  dimensions:  the  tortuosity  fractal  dimen¬ 
sion  Dt,  which  can  be  obtained  using  a  box-counting  method, 
and  the  pore  area  fractal  dimension  Df. 

The  pore  area  fractal  dimension  can  be  determined  by  sev¬ 
eral  methods,  which  can  be  roughly  sorted  into  three  types. 
Taking  linear  regression  of  the  experimental  data  to  gain  the 
fractal  dimension  is  the  first  method.  Such  experiments  include 
mercury  porosimetry  and  BET  absorption.  Processing  the  SEM 
images  or  the  images  from  other  imaging  technology  to  compute 
the  fractal  dimension  belongs  to  the  second  method.  The  most 
common  method  based  on  images  is  the  box-counting  method; 
furthermore,  a  difference  box-counting  method  and  a  Sierpinski- 
type  gasket  or  basket  method  can  also  be  used.  The  method  by 
analytical  solution  of  the  fractal  dimension  belongs  to  the  third 
approach. 

Using  a  box-counting  method  to  process  the  SEM  images  of 
two  samples:  TGP-H-060  carbon  paper  and  TGP-H-060  carbon 
paper  treated  with  PTFE,  we  get  the  tortuosity  dimension  of  two 
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Fig.  1.  A  comparison  of  permeability  from  the  present  fractal  model  and 
Kozeny-Carman  equation. 
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Fig.  2.  Effect  of  Df  and  Dj  on  permeability  for  carbon  paper. 


From  the  results,  it  can  be  seen  that  because  the  porosity  of 
carbon  paper  descends  and  the  pore  space  decreases  after  being 
treated  by  PTFE,  the  corresponding  permeability  decreases  with 
respect  to  that  of  carbon  paper  without  treatment. 

The  effects  of  the  pore  area  dimension  Df,  the  tortuosity 
dimension  Dx  and  porosity  0  on  the  permeability  are  plotted 
in  Figs.  2  and  3.  Fig.  2  shows  the  permeability  variation  with 
the  pore  area  dimension  Df  for  different  values  of  the  tortuosity 
dimension  Dj.  It  is  seen  that  the  permeability  increases  as  Df 
increases.  Increase  in  Df  corresponds  to  increase  in  the  number 
of  the  larger  pores,  and  decrease  in  the  population  of  the  smaller 
pores  within  a  representative  carbon  paper  volume.  The  result¬ 
ing  augmentation  in  the  available  area  for  gas  flow  in  turn  leads 
to  an  increase  in  the  permeability. 

An  interesting  fact  can  be  elucidated  from  Fig.  2  which  is 
that  for  any  tortuosity  dimension,  as  the  pore  area  dimension 
Df  approaches  its  largest  possible  value  of  2,  the  permeabil- 
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Fig.  3.  Effect  of  Dj  and  effective  porosity  on  permeability  for  carbon  paper. 


ity  approaches  maximum  possible  values  as  mentioned  above. 
This  theoretical  limit  on  Df  corresponds  to  the  situation  wherein 
the  carbon  paper  entirely  consists  of  pores  and  does  not  con¬ 
tain  any  solid  matrix.  Under  such  a  condition,  the  flow-through 
porous  medium  whose  porosity  approaches  unity  becomes  the 
flow  through  a  zone.  Moreover,  when  the  pore  area  dimension 
Df  approaches  its  lowest  possible  value  of  1 ,  the  permeability 
approaches  minimum  possible  values.  Its  theoretical  limit  on  Df 
corresponds  to  the  situation  wherein  the  carbon  paper  consists 
of  pores  of  sizes  approaching  zero,  which  no  longer  contributes 
to  the  gas  flow. 

The  variation  of  the  permeability  with  porosity  0  and  the  tor¬ 
tuosity  dimension  Dx  is  presented  in  Fig.  3.  The  permeability  is 
seen  to  decrease  with  increasing  the  tortuosity  dimension.  The 
decrease  in  permeability  with  the  increase  in  Dx  is  attributed  to 
the  increased  flow  resistance  due  to  the  highly  convoluted  cap¬ 
illary  pathways.  Furthermore,  the  permeability  increases  with 
rising  porosity  in  the  carbon  paper  owing  to  increase  in  the  avail¬ 
able  pore  volume  for  the  gas  to  permeate  through,  which  is  also 
observed  in  Fig.  2. 

The  effect  of  K 2  on  the  total  permeability  is  also  studied  with 
the  ratio  of  K2  to  the  total  permeability  K.  The  results  show  that 
the  ratios  lie  in  the  range  from  0.29  to  0.34%  for  carbon  paper 
and  from  0.33  to  0.38%  for  carbon  paper  treated  with  PTFE. 
The  values  of  K2IK  are  so  small  that  we  can  ignore  K2  in  the 
permeability  computation,  which  indicates  little  reaction  of  the 
Knudsen  effect  to  the  flow  through  carbon  paper,  so  the  modified 
items  of  Eqs.  (23),  (24)  and  (26)  need  not  be  taken  into  account 
in  permeation  through  carbon  paper.  In  addition,  the  ratios  of  K2 
to  the  total  permeability  K  of  the  number  two  sample  are  greater 
than  that  of  the  number  one  sample  for  the  mean  pore  size  of  the 
number  two  sample  is  smaller  than  that  of  number  one. 

5.  Conclusion 

In  this  paper,  the  relationship  between  the  microstructure  and 
the  permeability  of  the  GDL  in  a  PEMFC  was  given.  With 
the  microstructural  parameters  including  two  fractal  dimen¬ 
sions  of  porous  medium,  the  smallest  and  biggest  pore  diameter 
and  porosity  we  could  predict  the  corresponding  permeability 
through  fractal  models.  The  analytical  solution  permeability 
model  of  the  GDL  does  not  have  any  empirical  constants,  and 
every  parameter  in  this  model  has  a  specific  physical  meaning. 
So  the  permeability  model  gained  by  fractal  theory  is  obvi¬ 
ously  better  than  the  well-known  Konzeny-Carman  equation 
from  the  viewpoint  of  describing  a  penetrative  mechanism.  The 
prediction  results  show  that  the  present  fractal  permeability 
model  is  in  better  agreement  with  the  reference  value  than  the 
Kozeny-Carman  equation. 

The  effects  of  the  area  dimension  Df ,  the  tortuosity  dimension 
Dx  and  porosity  0  on  the  permeability  are  discussed  in  detail. 
And  we  can  draw  the  conclusion  that  permeability  increases  as 
Df  or  0  increases,  or  as  Dx  decreases. 

Therefore,  using  fractal  geometrical  theories  to  study  the  pore 
structures  of  the  GDL  in  a  PEMFC  provides  a  new  research 
method  for  the  gas  transport  mechanism  in  the  simulation  of 
PEMFCs.  It  makes  quantitative  analysis  of  the  relationship 
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between  the  micro-mechanism  and  the  macroscopical  material 
parameters  of  a  GDL  possible. 
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